In [1]:
import vmodel
import pygslib
In [2]:
parameters = {
'outfl' : None, # path to the output file or None (to use '_xxx_.out')
'nlag' : 40, # number of lags
'ivdir' : [[ 0.0,0.0,0.5 ],
[90.0,0.0,0.5]], # azm,dip,lag distance (array with shape [ndir,3])
'c0' : 0.2, # nugget effect
'it' : [1,1], # list of structure types (array with shape [nst])
'cc' : [0.4,0.4], # list of structure variances (array with shape [nst])
'ang1' : [0,0], # list of structure 1st rotation (array with shape [nst])
'ang2' : [0,0], # list of structure 2nd rotation (array with shape [nst])
'ang3' : [0,0], # list of structure 3rd rotation (array with shape [nst])
'a_hmax' : [10,10], # list of structure maximum horizontal ranges (array with shape [nst])
'a_hmin' : [5,5], # list of structure minimum horizontal ranges (array with shape [nst])
'a_vert' : [0,0]} # list of structure vertical ranges (array with shape [nst])
In [3]:
vg, fig, ax = vmodel.vmodel(parameters)
Parameters for VMODEL
*********************
START OF PARAMETERS:
_xxx_.out - file for variogram output
2 40 - number of directions and lags
0.0 0.0 0.5
90.0 0.0 0.5 - azm, dip, lag distance (array with shape [ndir,3])
2 0.2 - nst, nugget effect
1 0.4 0 0 0 - it,cc,ang1,ang2,ang3
10 5 0 - a_hmax, a_hmin, a_vert
1 0.4 0 0 0 - it,cc,ang1,ang2,ang3
10 5 0 - a_hmax, a_hmin, a_vert
VMODEL Version: 2.905
output file = _xxx_.out
ndir,nlag = 2 40
azm, dip, lag = 0.000000E+00 0.000000E+00 5.000000E-01
x,y,z offsets = 0.000000E+00 5.000000E-01 0.000000E+00
azm, dip, lag = 90.000000 0.000000E+00 5.000000E-01
x,y,z offsets = 5.000000E-01 6.083982E-09 0.000000E+00
nst, c0 = 2 2.000000E-01
it,cc,ang[1,2,3] = 1 4.000000E-01 0.000000E+00
0.000000E+00 0.000000E+00
a1 a2 a3 = 10.000000 5.000000 0.000000E+00
anis 1 2 = 5.000000E-01 0.000000E+00
it,cc,ang[1,2,3] = 1 4.000000E-01 0.000000E+00
0.000000E+00 0.000000E+00
a1 a2 a3 = 10.000000 5.000000 0.000000E+00
anis 1 2 = 5.000000E-01 0.000000E+00
VMODEL Version: 2.905 Finished
Stop - Program terminated.
In [4]:
vg
Out[4]:
average separation
var funct
number of directions
covariance
correlation
tail
head
type
cut
Variogram
Direction
Lag
0
0
1
0.0
NaN
2
1.00000
1.00000
None
None
None
None
1
0.0
0.20001
2
0.79999
0.79999
None
None
None
None
2
0.5
0.25995
2
0.74005
0.74005
None
None
None
None
3
1.0
0.31960
2
0.68040
0.68040
None
None
None
None
4
1.5
0.37865
2
0.62135
0.62135
None
None
None
None
5
2.0
0.43680
2
0.56320
0.56320
None
None
None
None
6
2.5
0.49375
2
0.50625
0.50625
None
None
None
None
7
3.0
0.54920
2
0.45080
0.45080
None
None
None
None
8
3.5
0.60285
2
0.39715
0.39715
None
None
None
None
9
4.0
0.65440
2
0.34560
0.34560
None
None
None
None
10
4.5
0.70355
2
0.29645
0.29645
None
None
None
None
11
5.0
0.75000
2
0.25000
0.25000
None
None
None
None
12
5.5
0.79345
2
0.20655
0.20655
None
None
None
None
13
6.0
0.83360
2
0.16640
0.16640
None
None
None
None
14
6.5
0.87015
2
0.12985
0.12985
None
None
None
None
15
7.0
0.90280
2
0.09720
0.09720
None
None
None
None
16
7.5
0.93125
2
0.06875
0.06875
None
None
None
None
17
8.0
0.95520
2
0.04480
0.04480
None
None
None
None
18
8.5
0.97435
2
0.02565
0.02565
None
None
None
None
19
9.0
0.98840
2
0.01160
0.01160
None
None
None
None
20
9.5
0.99705
2
0.00295
0.00295
None
None
None
None
21
10.0
1.00000
2
0.00000
0.00000
None
None
None
None
22
10.5
1.00000
2
0.00000
0.00000
None
None
None
None
23
11.0
1.00000
2
0.00000
0.00000
None
None
None
None
24
11.5
1.00000
2
0.00000
0.00000
None
None
None
None
25
12.0
1.00000
2
0.00000
0.00000
None
None
None
None
26
12.5
1.00000
2
0.00000
0.00000
None
None
None
None
27
13.0
1.00000
2
0.00000
0.00000
None
None
None
None
28
13.5
1.00000
2
0.00000
0.00000
None
None
None
None
29
14.0
1.00000
2
0.00000
0.00000
None
None
None
None
...
...
...
...
...
...
...
...
...
...
...
1
12
5.5
1.00000
2
0.00000
0.00000
None
None
None
None
13
6.0
1.00000
2
0.00000
0.00000
None
None
None
None
14
6.5
1.00000
2
0.00000
0.00000
None
None
None
None
15
7.0
1.00000
2
0.00000
0.00000
None
None
None
None
16
7.5
1.00000
2
0.00000
0.00000
None
None
None
None
17
8.0
1.00000
2
0.00000
0.00000
None
None
None
None
18
8.5
1.00000
2
0.00000
0.00000
None
None
None
None
19
9.0
1.00000
2
0.00000
0.00000
None
None
None
None
20
9.5
1.00000
2
0.00000
0.00000
None
None
None
None
21
10.0
1.00000
2
0.00000
0.00000
None
None
None
None
22
10.5
1.00000
2
0.00000
0.00000
None
None
None
None
23
11.0
1.00000
2
0.00000
0.00000
None
None
None
None
24
11.5
1.00000
2
0.00000
0.00000
None
None
None
None
25
12.0
1.00000
2
0.00000
0.00000
None
None
None
None
26
12.5
1.00000
2
0.00000
0.00000
None
None
None
None
27
13.0
1.00000
2
0.00000
0.00000
None
None
None
None
28
13.5
1.00000
2
0.00000
0.00000
None
None
None
None
29
14.0
1.00000
2
0.00000
0.00000
None
None
None
None
30
14.5
1.00000
2
0.00000
0.00000
None
None
None
None
31
15.0
1.00000
2
0.00000
0.00000
None
None
None
None
32
15.5
1.00000
2
0.00000
0.00000
None
None
None
None
33
16.0
1.00000
2
0.00000
0.00000
None
None
None
None
34
16.5
1.00000
2
0.00000
0.00000
None
None
None
None
35
17.0
1.00000
2
0.00000
0.00000
None
None
None
None
36
17.5
1.00000
2
0.00000
0.00000
None
None
None
None
37
18.0
1.00000
2
0.00000
0.00000
None
None
None
None
38
18.5
1.00000
2
0.00000
0.00000
None
None
None
None
39
19.0
1.00000
2
0.00000
0.00000
None
None
None
None
40
19.5
1.00000
2
0.00000
0.00000
None
None
None
None
41
20.0
1.00000
2
0.00000
0.00000
None
None
None
None
84 rows × 9 columns
In [ ]:
Content source: opengeostat/pygslib
Similar notebooks: